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ABSTRACT 

We develop simple analytic representations of the polarized quark and gluon 
distributions in the nucleon at low Q^ which incorporate general constraints ob- 
tained from the requirements of color coherence of gluon couplings at a; ~ and 
the helicity retention properties of perturbative QCD couplings at a; ~ 1. The un- 
polarized predictions are similar to the D'^ distributions given by Martin, Roberts, 
and Stirling. The predictions for the quark helicity distributions are compared 
with polarized structure functions measured by the E142 experiment at SLAC and 
the SMC experiment at CERN. 



1. INTRODUCTION 

Measurements of polarization correlations in high momentum transfer reactions 
can provide highly sensitive tests of the underlying structure and dynamics of 
hadrons. The most direct information on the light-cone momentum distributions 
of helicity-aligned and helicity-anti-aligned quarks in nucleons is obtained from 
deep inelastic scattering of polarized leptons on polarized targets. Recent fixed- 
target measurements, including the CERN SMC muon-deuteron experiment , the 
electron- iJe"^ and electron-proton experiments E142 and E143 at SLAC , and the 
SMC muon-proton experiment are now providing important new constraints on 
the proton and neutron helicity-dependent structure functions. 

Although the Q^-evolution of structure functions is well-predicted by pertur- 
bative QCD, the initial shape of these distributions reflects the non-perturbative 
quark and gluon dynamics of the bound-state solutions of QCD. Nevertheless, it is 
possible to predict some aspects of the shape of the input nucleon structure func- 
tions from perturbative arguments alone. In this paper, we will develop simple 
analytic representations of the quark and gluon helicity distributions which incor- 
porate general constraints obtained from the color coherence of the gluon couplings 
at X ~ 0, and the helicity structure of perturbative QCD couplings at a; ~ 1. Since 
we work at the bound-state scale, we can directly impose global sum rules and 
symmetries such as the axial coupling constraint A-u = Arf = qa^ from neutron 
beta decay. The parameterizations we use have the minimal number of parameters 
needed to satisfy all of the constraints. The predicted forms for the quark and 
gluon helicity distributions Ag(x) = q^{x) — <i~{x) and AG(a;) = G^{x) — G^{x) 
should provide useful guides to the expected shapes of the polarized structure func- 
tions and an understanding of how the helicity content of the nucleon is distributed 
among its constituents. Eventually these fundamental distributions should be com- 
putable using non-perturbative methods such as lattice gauge theory or light-cone 
Hamiltonian diagonalization. 

The structure functions discussed in this paper are meant to reflect the intrinsic 



bound state structure of the nucleons, and thus they strictly apply only at low 
resolution scales Q^ < Qq where QCD evolution can be neglected. They can be 
used at large Q^ as the input distributions for perturbative QCD evolution, as in 
the analyses of Ref . [4] . However, as we discuss below, the actual implementation 
the evolution program must take into account the fact that for a; ~ 1 the bound 
state quark which is struck in deep lepton inelastic scattering is far off shell, thus 
suppressing its gluon radiation. 

The polarized quark and gluon distributions Gq/jj{x,X,Q) and Gg/ff{x,X,Q) 
of a hadron are most simply represented as probability distributions determined by 
the light-cone wavefunctions ipn{xi,k_i_i, \i), where XlILi -^i ~ -'-' ^^^ SiLi ^J-i ~ 
0^. The square of the invariant mass of an n-particle Fock State configuration in 
the wavefunction is Ai'^ = X]iLi(^i_i"'""^?)/-^*- Thus the kinematical regime where 
one quark has nearly all of the light-cone momentum a; ~ 1, and the remaining 
constituents have Xi ~ 0, represents a very far off-shell configuration of a bound 
state wavefunction. In the limit x — *> 1, the Feynman virtuality of the struck parton 
in a bound state becomes far off-shell and space-like: k'^ — m^ = x{Mfj — A4'^) -^ 
—fi^/{l — x), where fx is the invariant mass of the system of stopped constituents. If 
one assumes that the bound state wavefunction of the hadron is dominated by the 
lowest invariant mass partonic states, then the constituents can attain far off-shell 
configurations only by exchanging hard gluons; thus the leading behavior at large 
virtuality can be computed simply by iterating the gluon exchange interaction 
kernel ^' ' . This conforms to the usual ansatz of perturbative QCD that hard 
perturbative contributions dominate amplitudes involving high momentum transfer 
compared to the contributions arising from non-perturbative sources. 

Thus, because of asymptotic freedom, the leading order contributions in as{k'^) 
to the quark and gluon distributions at a; ^ 1 can be computed in perturbative 
QCD from minimally connected tree graphs. For example, in the case of the 
nucleon structure functions, the dominant amplitude is derived from graphs where 
the three valence quarks exchange two hard gluons. The tree amplitude is then 
convoluted with the nucleon distribution amplitude 0(a;j,A;^) which is obtained 



by integrating the valence three-quark nucleon wavefunction ilj^{xi,k±i,Xi), over 
transverse momenta up to the scale k'j^ . The dk±^d(j) azimuthal loop integrations 
project out only the L^ = component of the three-quark nucleon wavefunction. 
Thus, in amplitudes controlled by the short distance structure of the hadron's 
valence wavefunction, orbital angular momentum can be ignored, and the valence 
quark helicities sum to the hadron helicity. 

The limiting power-law behavior at a; ^ 1 of the helicity-dependent distribu- 
tions derived from the minimally-connected graphs is 



Gq/H ~ (1 - a;) 



where 

p = 2n - 1 + 2^S^. 

Here n is the minimal number of spectator quark lines, and AS'^ = \Sl — S^\ = 
or 1 for parallel or anti-parallel quark and proton helicities, respectively ' . This 
counting rule reflects the fact that the valence Fock states with the minimum num- 
ber of constituents give the leading contribution to structure functions when one 
quark carries nearly all of the light-cone momentum; just on phase-space grounds 
alone, Fock states with a higher number of partons must give structure functions 
which fall off faster at a; ^ 1. The helicity dependence of the counting rule also re- 
flects the helicity retention properties of the gauge couplings: a quark with a large 
momentum fraction of the hadron also tends to carry its helicity. The antiparallel 
helicity quark is suppressed by a relative factor (1 — x)^. Similarly, in the case of 
a splitting function such as q ^ qg ot g ^ qq, the sign of the helicity of the par- 
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ent parton is transferred to the constituent with the largest momentum fraction . 
The counting rule for valence quarks can be combined with the splitting functions 
to predict the x ^ 1 behavior of gluon and non-valence quark distributions. In 
particular, the gluon distribution of non-exotic hadrons must fall by at least one 
power faster than the respective quark distributions. 



The counting rules for the end-point-behavior of quark and gluon hehcity dis- 
tributions can also be derived from duality, i.e., continuity between the physics 
of exclusive and inclusive channels at fixed invariant mass ' . As shown by Drell 
and Yan , a quark structure function Gg/jj ~ (1 — a;)^""^^ at a; ^ 1 if the cor- 
responding form factor F(Q ) ~ (l/Q )"" at large Q . Recent measurements of 
elastic electron-proton scattering at SLAC are compatible with the perturbative 
QCD predictions for both the helicity-conserving Fi{Q'^) and helicity-changing 
F2{Q^) form factors: Q^Fi{Q'^) and Q^F2{Q'^) become approximately constant at 
large Q^. The power-law fall-off of the form factors corresponds to the helicity- 
parallel and helicity-antiparallel quark distributions behaving at a; ^ 1 as (1 ~ x)^ 
and (1 — x)^, respectively, in agreement with the counting rules. The leading ex- 
ponent for quark distributions is odd in the case of baryons and even for mesons 
in agreement with the Gribov-Lipatov crossing rule 

The counting rule predictions for the quark and gluon distributions are relevant 
at low momentum transfer scales Qq ~ ^qcd in which the controlling physics is 
that of the hadronic bound state rather than the radiative corrections associated 
with structure function evolution. At the hadronic scale one can normalize the 
non-singlet quark helicity content of the proton and neutron using the constraint 
from P decay : 

Au- Ad= — = 1.2573 ± 0.0028 . 
9V 

where Aqi{x) = (il{x)—q^{x) with i = -u, rf, s is the difference of the helicity-aligned 

and helicity- anti- aligned quark distributions in the proton, and Aqi = Jq dxAqi{x) 

is the integrated moment. [In the standard notation q^{x,Q) = Gq/p{x,Xq = 

Ap, Q) + Gq/p{x, Xq = Xp, Q) so that both quark and anti-quark contributions are 

included.] In addition, if one assumes SU{3) flavor symmetry, hyperon decay also 

implies a polarized strange quark component in the proton wavefunction 

^" + ^^-^^-^ = 0.39 . 

Thus only one normalization is left undetermined. 



The presence of polarized gluons in the nucleon wavefunction imphes that 
polarized strange quarks contribute to the nucleon helicity-dependent structure 
functions at some level. There is also evidence from neutrino-proton elastic scat- 
tering that the proton has a significant polarized strange quark content ' . Our 
parameterization of the polarized strange quark distributions with a significant he- 
licity fraction As ~ 0.10 and a small momentum fraction (xg) — 0.035. The shape 
represents the sum of contributions from radiatively-generated s + s quarks as well 
as intrinsic strange quarks intrinsic to the nucleon bound states. 

The helicity-dependent structure function gi{x,Q^) measured in deep inelastic 
polarized-lepton polarized-proton scattering can be identified in the Bjorken scaling 
region with the quark helicity asymmetry: 



The first moment of the proton-neutron difference has zero anomalous dimension 
and satisfies the Bjorken sum rule: 

1 

J dx[g{{x, Q2) _ ,^(a:, q2)] = ^ ^ (l - ^ + ■■■), 


where the last factor represents the radiative corrections from hard gluon interac- 
tions in the electron-quark scattering process . Thus the QCD radiative correc- 
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tions to the helicity-dependent structure functions can modify the shape of the 
distributions, within the global constraint of the Bjorken sum rule. 

At high Q^, the radiation from the struck quark line increases the effective 
power law fall-off (1 — x)^ of structure functions relative to the underlying quark 
distributions: Ap = {4CF/(3i)log{logQ^/A^)/{\ogQl/A^) where Cp = 4/3 and 
f3i = 11 — (2/3)n/. The counting rule predictions for the power p thus provide a 
lower bound for the effective exponent of quark structure functions at high Q^ > Qq. 
However, in the end-point region a; ~ 1, the struck quark is far off-shell and the 



radiation is quenched since one cannot evolve Q below Q\':^kh = — (/i /(I — x)\ 
the Feynman virtuality of the struck parton . Furthermore, the integral of the g\ 
structure function is only affected by QCD radiative corrections of order as((5^)/7r. 

Thus PQCD can give useful predictions for the power law fall-off of helicity- 
aligned and anti- aligned structure functions at a; ~ 1. Higher order contributions 
involving additional hard gluon exchange are suppressed by powers of q;s(A;^). 
Further iterations of the interaction kernel will give factors of fractional powers 
of log(l — x) analogous to the anomalous dimensions log'*'" Q^ which appear in 

fl2l 

the PQCD treatment of form factors at large momentum transfer . This is in 
contrast to super-renormalizable theories such as QCD (1+1) where the power-law 
behavior in the endpoint region is modified by all-order contributions 

The fact that one has a definite prediction for the a; ~ 1 behavior of leading 
twist structure functions is a powerful tool in QCD phenomenology, since any 
contribution that does not decrease sufficiently fast at large x is most likely due to 
coherent multi-quark correlations. As discussed in Ref. [21], such contributions are 
higher twist, but they arise naturally in QCD and are significant at fixed (1— a;)(5^. 
Such coherent contributions are in fact needed in order to explain the anomalous 
change in polarization seen in pion-induced continuum lepton-pair and hadronic 
J/-?/' production experiments at high xp 

At large x the perturbative QCD analysis predicts "helicity retention" - ie., 
the helicity of a valence quark with x ~ 1 will matches that of the parent nucleon. 
This result is in agreement with the original prediction of Farrar and Jackson that 
the helicity asymmetry I\q(x) approaches 1 at a; ^ 1. We also predict, in agreement 
with Ref. [6], that the ratio of unpolarized neutron to proton structure functions 
approaches the value 3/7 for x — > 1. 

In the following sections we will analyze the shape of the polarized gluon and 
quark distributions in the proton. First we will study the behavior of the gluon 
asymmetry I\G{x)IG{x) (polarized over unpolarized distributions) at small val- 
ues of x, where it turns out to be proportional to x with a coefficient approxi- 



mately independent of the details of the bound-state wavef unction. We then write 
down a simple model for the gluon distributions which incorporates the count- 
ing rule constraints at a; — > 1. The same is done for the up, down and strange 
quark distributions. The extrinsic and intrinsic strange quark distributions are 
also discussed, paying special attention to the inclusive-exclusive connection with 
the strange quark contribution to the proton form factors. 

2. HELICITY-DEPENDENT GLUON DISTRIBUTIONS 

The angular momentum of a fast-moving proton has three sources, the angular 
momentum carried by the quarks, the angular momentum carried by the gluons, 
and the orbital angular momentum carried by any of the constituents. Angular 
momentum conservation for J^ at a fixed light-cone time implies the sum rule 

i(AM + Arf + As) + AG' + <L^ >= -. (2.1) 

Here AG = Jq dxAG{x) is the helicity carried by the gluons, where AG{x) is the 
difference between the helicity-aligned and anti- aligned gluon distributions G^{x) 
and G^{x); the unpolarized gluon distribution G{x) is the sum of these two func- 
tions, G{x) = G'^lx) + G^{x). The corresponding definitions for the quark dis- 
tributions Aq{x) = q^{x) — q^{x) and q{x) = q^{x) + q^{x) with q = u,d,s. By 
definition, the antiquark contributions are included in Aq{x) and q{x). As empha- 
sized by Ma , the helicity distributions measured on the light-cone are related 
by a Wigner rotation (Melosh transformation) to the ordinary spins Sf of the 
quarks in an equal-time rest-frame wavefunction description. Thus, due to the 
non-colinearity of the quarks, one cannot expect that the quark helicities will sum 
simply to the proton spin. 

In this paper we shall present model forms for the gluon distribution functions 
AG{x) and G{x) for nucleons which incorporate the known large-x counting-rule 



constraints: 

G+{x)^C{l-x)^ (x^l), (2.2) 

G-{x)^C{l-xf (x^l); (2.3) 

We will also incorporate a constraint on the behavior of the gluon asymmetry ratio 
AG{x)/G{x) for small x : 

fAG{x)\ X l\ 

(^) /proton^ 3 \l/ 



\J^\ ^^{-.) (--0). (2.4) 



This last theoretical constraint will be demonstrated below. Here (l/y) stands for 
the first inverse moment of the quark light-cone momentum fraction distribution 
in the proton lowest Fock state. For this state we expect (l/y) — 3. 

A simple form for baryon gluon distributions, which incorporates the limiting 
behaviors presented above, is 

AG{x)=- [l-{l-xf] (l-x)^ 

N (2.5) 

G(x)=- [l + n-x?] (1-x)^. 

X 

In this model the momentum fraction carried by the gluons in the proton is (xg) = 
L dxxG{x) = ^N, and the helicity carried by the gluons is AG = Jq dxAG{x) = 
(11/30)A^. Taking the momentum fraction (xg) to be 1/2, we predict AG = 0.54. 

Such large values for the gluon momentum fraction are inconsistent with the 
assumption that the proton has a dominant three-quark Fock state probability; 
a self-consistent approach thus requires taking into account gluon radiation from 
the full quark and gluon light-cone Fock basis of the nucleon. Our main emphasis 
here is to predict the characteristic shapes of the polarized quark and gluon dis- 
tributions. The large x regime is clearly dominated by the lowest particle-number 
Fock states. We thus expect that the qualitative features of the model to survive 
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in a more rigorous approach; in particular, it is apparent from the structure of the 
model, that the gluon helicity fraction will be of the same order of magnitude as 
the gluon momentum fraction. 

The prediction that AG ^ 0.5 is phenomenologically interesting. If one also 
accepts the experimental suggestion from EMC that the quark helicity sum Au+Ad 
is small, then this implies that gluons could carry a significant fraction of the proton 
helicity J^ = 1/2 of the same size as the momentum fraction carried by the gluons. 
However, one also expects significant orbital angular momentum Lz which arises, 
for example, from the finite transverse momentum associated with the q -^ qg 
gluon emission matrix element. 

We now proceed to prove equation (2.4) for the low-x behavior of the asym- 
metry AG{x)/G{x). In this region the quarks in the hadron radiate coherently, 
and we must consider interference between amplitudes in which gluons are emitted 
from different quark lines. An analysis of this type was first presented in Ref. [25] , 
and in this note we extend and correct some of the results of that paper. 

As an example, we first analyze the helicity content of positronium, where we 
can ignore internal transverse momenta and non-collinearity. Consider the ortho- 
positronium two-fermion Jz = 1 Fock state in which the particles have helicities 
+ +. Following the calculation of Ref. [25], we obtain 

f^^'^''^^ r^x(-\c^2x (x^O). (2.6) 



V c!{x) J ortho{j,^+i) \y 

In the case of para-positronium (and also for J^ = ortho-positronium), in which 

we start with a Fock state with helicities H , the result is AG{x) = 0. This is 

because for every diagram in G'+(a;) there is a corresponding diagram in G~{x), 
but with the helicities of all the particles reversed. 

We now apply a similar analysis to the gluon distribution in the nucleon. We 
start with a three-quark Fock state in which the quarks have helicities + + + as 
would be appropriate for the helicity content of an isobar state A with Jz = 3/2. 
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Then the resuh found in Ref. [25], i.e. 



'^(^) /A(J,=3/2) \^/ 



follows. 



In the nucleon case, however, we start with a three-quark Fock state with 
helicities + + — . Thus clearly there is a cancellation between the squared terms 
in which the gluon is emitted from one of the positive helicity quarks versus the 
contributions in which the gluon is emitted by a negative helicity quark. The 
interference terms work similarly, ensuring a finite result for both G{x) and AG{x) 
at zero kj_, just as in the case of photon distributions in positronium. Then the 
positive helicity quarks have a dominant G'+(a;) and contribute positively to AG{x); 
similarly, the negative helicity quarks contribute negatively to AG{x). To see this 
more clearly, consider the photon emitted by a single electron with J^ = +1/2. 
Then G^, (x) = 1/x and G^,{x) = (1 — x)^ jx. Thus AG{x) / G{x) = a; at a; -^ 
with unit coefficient in this case. The sign reverses for an electron with Jz = —1/2. 

The generated gluon asymmetry distribution in the nucleon at low x is then 
given by equation (2.4). The extra factor of 1/3 is due to the fact that all the 
quarks contribute positively to G{x), but they give contributions proportional to 
the sign of their helicity in AG{x). The main assumption setting the value of 
the gluon asymmetry at a; ^ is the estimated value of the inverse moment 
(l/y). For realistic wavefunctions this expectation value may receive very large 
(possibly divergent) contributions from near y = 0. However, one must be careful 
at this point because in deriving Eq. (2.4) we assumed that x <^ y. In order 
to be consistent with this assumption we will perturb around a constituent quark 
wavefunction which is strongly peaked around y = (y) = 1/3. We have furthermore 
assumed for simplicity that (y) is the same for all valence quarks, although this is 
inconsistent with results from QCD sum rules . (One could improve the estimate 
for {l/y) by allowing for different momentum fractions for the helicity-up and 
helicity-down quarks. This would evidently reduce AG, since it is known that {y) 
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is larger for helicity parallel quarks. Furthermore, in QCD we expect that higher 
Fock states will contribute to reduce the value of {y) away from 1/3, which would 
be the expected value if only the three-quark valence Fock state was present.) 

3. HELICITY-DEPENDENT QUARK DISTRIBUTIONS 

In this section we shall construct a simple polynomial model for the helicity- 
dependent quark distributions in the proton and neutron. 

As we have discussed in the previous sections, at a; ~ 1 PQCD predicts that 
the helicity-parallel quark distribution q^{x) is enhanced relative to the helicity- 
antiparallel quark distribution q~{x) by two powers of (1 — x). The property of 
helicity retention at large a; is a direct consequence of the gauge theory couplings 
between quarks and gluons. For the valence quarks in a nucleon the counting rules 
predict 

g+(x) ~ (l-a;)3 (x ^ 1), (3.1) 

and 

q~{x)^{l-xf (x^l). (3.2) 

The case of the non-valence strange quarks is somewhat more complex and will be 
discussed in detail in the next section. The result is 

s^{x)^{l-xf (x^l), (3.3) 

s-{x)^{l-xf {x^l). (3.4) 

For X ~ the helicity correlation disappears since the constituent has infinite 
rapidity Ay = logx relative to the nucleon's rapidity. 
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The strange quark distribution in a nucleon can arise from both intrinsic and 
extrinsic contributions. The intrinsic contribution is associated with the muhipar- 
ticle Fock state decomposition of the hadronic wavefunction, and it is essentially 
of non-perturbative origin. This is in contrast to the extrinsic component, which 
arises from ss pair production from a gluon emitted by a valence quark, and is 
associated with the self-field of a single quark in the proton. From evolution and 
gluon splitting, the extrinsic strange contributions are known to behave as 

st{x) r^ {l-xf (x^l) (3.5) 

s~{x)r^{l-xy (x^l). (3.6) 

The Drell-Yan inclusive-exclusive connection relates the high Q^ behavior of 
the hadronic form factors to the large x limit of the quark distribution functions; 
i.e. 

i^(Q^)^^^(^ ^ G,/,^-^(l-.)2-^+2A..^ (3,7) 

where ASz = or 1 for parallel or antiparallel quark and proton helicities, re- 
spectively. If we naively apply this prescription to the extrinsic strange quark 
component, we would predict that the strange quark contribution to the electro- 
magnetic proton form factor should fall as l/Q^, since in this case n = 3. But 
a direct calculation of the strange quark contribution to either the axial or vec- 
tor form factor of the nucleon gives only a nominal 1/Q^ behavior, which is the 
same power-law fall-off as the valence quark contribution. In the leading order 
calculations the loop integrals connecting a hard ss loop to a valence quark all 
have momenta i = 0{Q), thus producing radiative corrections of order ag{Q), to 
the exclusive amplitude with N = 2 (axial) or A^ = 3 (vector), rather than extra 
powers of l/Q'^ . The solution to this apparent contradiction is that we should 
apply the inclusive-exclusive connection for the strange quark contributions to a 
transition form factor connecting an initial state with three quarks (uud) to a final 
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state in which a strange pair has been created {uudss), as in the transition form 
factor p — » AK, at fixed final state mass. Since the internal hard-scattering matrix 
element Tjj for (uud) + 7* -^ sudus has three off-shell fermion legs, this transition 
form factor falls off as (l/Q ) , and it correctly satisfies the inclusive-exclusive 
connection (n = 3). 

One can also consider the case where Q^ and the final state mass are both large, 
but there is a -ft' and A in the final state. This again corresponds to a ~ (1 — a;)^ 
structure function. In the case of the transition p -^ p(j), there is a color mismatch 
in Tfj at lowest order. Thus this amplitude should be suppressed (Zweig rule) by 
an extra power of a<,.((5^). Of course all of this holds for the analogous charm 
systems as well. 

The intrinsic strange components are associated with Fock states having at 
least five particles; the distributions thus have the behavior 

sf{x) r^ {l-xf (x^l) (3.8) 

s^{x)r^{l-x)^ (x^l), (3.9) 

which corresponds to n = 4 in the spectator quark counting rules. It also satisfies 
the inclusive-exclusive connection, since the intrinsic contribution to the form factor 
falls as (l/g2)4. 

For the complete parameterization we shall adopt the canonical forms: 



«''(^) = ^ [^-(1 - ^)^ + ^-(1 - ^)'] ' 



d+{x) = ^[Adil-xf + Bdil-x)^], 



U [X] 



l-[Cu{l-xf + Du{l-xf], 
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^[Cs{l-xf + Ds{l-xf], 



S [X 



where we require 



Aq + Bq = Cq + Dq 



to ensure the convergence of the hehcity-dependent sum rules. Thus in our model, 
the Regge behavior of the asymmetry Ag(a;) ~ x~"^ is automatically one unit less 
than the unpolarized intercept: a^ = a — 1. Isospin symmetry at low x (Pomeron 
dominance) also requires 

We emphasize that these distributions include both the quark and antiquark con- 
tributions. 

Our parameterization of the helicity-dependent quark distributions is close in 
spirit to the parameterization Dq for the unpolarized quark and gluon distributions 
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given by Martin, Roberts and Stirling . The MRS parameterizations are a good 
match to our unpolarized forms q{x) = q^{x)+q^{x) since the MRS forms combine 
counting-rule constraints with a good fit to a wide range of perturbative QCD 
phenomenology. We find that choosing the effective QCD Pomeron intercept a = 
1.12 allows good match to the unpolarized quark distributions given by the MRS 
parameterization Dq at Q^ = 4 GeV^ over the range 0.001 < x < 1. It also predicts 
an increasing structure function F2{x,Q ) for x < 10^ , as suggested the recent 

[281 ^^^ 

data from HERA . Thus we predict aR = 0.12 for the helicity-changing Reggeon 
intercept. The momentum fraction carried by the quarks (and antiquarks), {xq) = 
Jq dxxq{x), where q{x) = q^{x) + q^{x), is assumed to be ~ 0.5. 
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A combined analysis '' of the SLAC and EMC ' polarized electron-proton 
data provides the constraint: 

/ dxg^{x) = 0.112 ± 0.009 ± 0.019 . 

If one uses the central value together with the constraints from nucleon and hy- 
peron decay and includes radiative corrections of 0{{as/7i)^) then one obtains the 
following values for the proton helicity carried by the different quarks ^ : 

Am = 0.83 ±0.03, Arf = -0.43 ± 0.03, As = -0.10 ± 0.03 , (3.10) 

at the renormalization scale Q^ = 10 GeV . Since these values for the Ag are 
obtained after removing the deep inelastic radiative corrections, we can use them 
as the initial phenomenological inputs for the proton; the neutron distributions 
then follow from isospin symmetry. The small value for the total quark helicity 
AS = Am + Ad + As = 0.31 ± 0.07 is consistent with large Nc predictions in 
QCD , and it is about half of the value AS ~ 0.75 expected in the relativistic 
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three-quark constituent model for the nucleon without dynamical gluons . As 
we shall show below, the gluon helicity fraction Ag scales closely with the gluon 
momentum fraction (xg) . 

The u{x) and d{x) parameterizations have eight parameters which we will fix 
using the following eight conditions: three conditions arise from the requirement 
that the sum rules converge at a; -^ 0; two conditions come from the values of 
Am and Ad; one condition follows by imposing the SU{6) large x relation Au = 
5Ad; one condition is obtained from the empirical value of the Gottfried sum 
Sg = J[u{x) — d{x)]/3 = 0.235 ; the final condition is obtained from the sum of 

[271 

momentum fractions carried by the quark and antiquark Xu + x^ = 0.521 . It is 
straightforward to find parameters for the polynomial forms which are consistent 
with the above inputs. 

A^ = 3.784, Ad = 0.757, (3.11) 
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Bu = -3.672, Bd = -0.645, (3.12) 

Cu = 2.004, Cd = 3.230, (3.13) 

L>^ = -1.892, Dd = -3.118,. (3.14) 

With this set of parameters, the respective quark momentum fractions are: 

{xu) = 0.331, (xd) = 0.190 (3.15) 

The predicted distributions xu{x) , xd{x) , Ad{x) , and Au{x) are shown in Figs. 
1(a) and 1(b). In each case both the quark and antiquark contributions are in- 
cluded. The simple polynomial forms represent a simple parameterization of the 
non-perturbative polarized and unpolarized quark distributions which satisfy the 
known theoretical constraints at large and small x and the empirical sum rules. 
We also show a comparison of the unpolarized distributions with the MRS Dq 
phenomenological parameterizations. The agreement is quite reasonable. The 
differences in the shapes of the distributions can be attributed to the effects of 
perturbative QCD evolution. 

Notice that Ad{x) = d^{x) — d^{x) is positive at large x (which follows from 
Alt = ^Ad), and negative at small to moderate values of x. We thus predict that 
Ad{x) will change sign and go through zero at some physical value for x. With 
the above parameterization the zero of Ad{x) occurs at a; = 0.489. 



In the case of the strange quark plus strange antiquark distributions, we have 
four parameters and three conditions: one from the convergence of sum rules; one 
from the value of As; and one from the momentum fraction carried by strange 
plus anti-strange quarks Xs = 0.035 . This leaves us with one unknown, which 
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we choose to be Cg- The three constraints give the solution set: 

A, = -0.6980 + 0.9877C,., Bs = 0.8534 - l.lUlCs, Ds = 0.1551 - 1.1294(7^ 

(3.16) 
Because of the probabihstic interpretation of parton distribution functions, s^(x) 
and s~{x) both must be non-negative functions of x, which imphes the rather 
stringent bounds 

0.7067 < Cs < 1.2013. 

Within these bounds, gi{x) is practically independent of Cs] to be definite, we chose 
Cs = 1. (We could have taken any other value consistent with the inequalities ^ .) 
We compare our simple parameterization to the MRS Dq parameterization in Fig. 
1(c). The MRS distribution which gives an approximate realization of the data 
rises faster at low x than our model. This could be attributed to the need to impose 
higher Pomeron intercept, or the effects of QCD evolution. 

We can also find parameterizations for the polarized gluon distributions which 
are consistent with the x ^ and a; — > 1 helicity constraints, as well as the MRS 
unpolarized gluon distribution: 

G+ix) = -L [A,il - xf + B,{1- xf] , (3.17) 

G-{x) = -^ [Ag{l - xf + Bg{l - xf] , (3.18) 

This form automatically incorporates the coherence constraint, Eq. (2.4). We shall 
assume that Ug = a = 1.12 so that the pomeron intercept is identical for quark 
and gluon distributions. The parameter set Ag = 2 and Bg = —1.25 gives an unpo- 
larized gluon distribution G{x) = G'^{x) + G^{x) similar to the phenomenological 
Dq gluon distribution given by MRS. (See Fig. 2.) The momentum carried by the 
gluons in the nucleon using the above simple form is (xg) = 0.42. (The gluon and 
light quark and antiquark distributions then almost saturate the momentum sum 
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rule.) The gluon helicity content for the above parameterization is AG = 0.45. As 
shown in the figure, the shape of the polarized distribution AG{x) given by the 
above parameterization is almost identical to xG{x). 

Alternatively, if we take ag = 1, then the parameter set Ag = 0.2381 and 
Bg = 1.1739. give again the same values (x — g) = 0.42 and AG = 0.45 as above. 
In this case, the resulting shape unpolarized distribution G{x) = G'^{x) + G^{x) is 
indistinguishable from the phenomeno logical Dq gluon distribution given by MRS. 

Although there is some experimental information about the unpolarized gluon 
distribution, this is not the case for the polarized gluon distribution. It is important 
to test these distributions directly, for example in processes such as J/ip production 
in polarized e — p and p — p collisions 

4. POLARIZED STRUCTURE FUNCTIONS 

In this section we will use the polynomial model forms for Aq{x) and q{x) to 
compute the polarized helicity structure functions of nucleons: 

g^ix) = \ {^Au{x) + ^-M{x) + \as{x)^ (4.1) 

and 

gT{x) = \ (^-M{x) + \au{x) + \as{x)^ , (4.2) 

and compare the results to the recent experiments. (Note that Aq{x) refers to the 
combined asymmetries from both quarks and antiquarks in the proton.) A precise 
prediction requires the inclusion of QCD evolution. Here we will, as in Ref. [15], 
simply include the normalization factor Nqcjj = 1 — (as/yr) ^ .92 arising from 
QGD radiative corrections. The Bjorken sum rule for the difference of proton and 
neutron quark helicities is automatically satisfied. The Ellis- Jaffe sum rule for the 
nucleon quark helicity is violated by the model due to the presence of the strange 
quark contributions As. 
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We have emphasized that the dynamics of QCD imphes hehcity retention: the 
quark with x close to 1 has the same hehcity as the proton. Thus all of the structure 
functions asymmetries become maximal at a: -^ 1, and the ratio of unpolarized 
proton and neutron structure functions can be predicted. 

According to the standard SU{6) flavor and helicity symmetry, the probabilities 
to find u and d quarks of different helicities in the proton's three-quark wavefunc- 
tion are: P{u+) = 5/9, P{d+) = 1/9, P{u-) = 1/9, P(rf-) = 2/9'"'. Thus the 
usual expectation from SU(6) symmetry is F2{n)/F2{p) = 2/3 for all x. As Farrar 
and Jackson pointed out , this naive SU{6) result cannot apply to the local he- 
licity distributions since the helicity aligned and helicity anti-aligned distributions 
have different momentum distributions. At large x -u" and d^ can be neglected 
relative to u'^ and (i+, and thus SU{6) is broken to S'f/(3)+ x SU{3)^ . Our model 
retains the SU{6) ratio P{u^)) : P{d^)) = Au '■ A^ = 5 : 1, at large x so that we 
predict F2{n)/ F2{p) -^ 3/7 as a; ^ 1. The physical picture that emerges is that the 
struck quark carries all the helicity of the nucleon, and the spectators have 5*^ = 0, 
although their total helicity is a combination of and 1. This wavefunction is just 
a piece of the full SU{Q) wavefunction, but since it is the piece that contains the 
u^ and d'^ , and this part remains unchanged, the ratio P{u^)/ P{d'^) is still 5/1. 

Notice that the only empirical input into our model is the integrated values 
of the various flavors obtained from the proton data. The shape of the polarized 
distributions is essentially determined by the perturbative QCD constraints. The 
agreement with the shape of the SLAC and EMC experimental data for the proton 
is rather good (see Fig. 3(a)) and could be further improved by taking into account 
PQCD evolution. 

We can also compare our model with the polarized neutron structure function 
extracted by the E142 from its polarized electron polarized He measurements. 
(See Fig. 3(b).) For the neutron we predict two new effects which are not present 
in the proton. First g^^ tends to fall faster than g'^ for large x. This is because 
as in the Carlitz-Kaur and Farrar- Jackson models, the helicity aligned up- 
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quark dominates the proton distribution and the hehcity down quark dominates 
the neutron structure function at large x. A related effect is that gf^ix) changes 
sign as a function of x. This is due to the fact that except for large x (where 
the helicity aligned down quark dominates) gf^ is dominated by the anti-aligned 
up quark distribution. Since Jq dxAuri{x) = Jq dxAdp{x) < "^ it is clear that 
gf^i^x) must be negative at small x. 

A comparison of our model with the recent SMC data for the polarized deuteron 
structure function gi{x) is shown in Fig. 4. The shape of the data appears to be 
consistent with our predictions, except possibly at the largest x point where the 
model predicts too little asymmetry. To make this prediction we have, as in Ref. 
[1], assumed that the deuteron structure function is half of the sum of the neutron 
and proton structure functions and included the D-state depolarization factor with 
D-state probability 0.058. The model then predicts the normalization 



I dxgfix) = \ \ dx{g\{x) ^ glix)) 



5 1 
— (Am + Arf) H As 



a^\ ( 3 

- -CUD = 0.038 



TT 



compared to the SMC result 

dxgfix) = 0.023 ± 0.020(stat.) ± 0.015(syst.) . 



The distributions presented in this paper have applicability to any PQCD 
leading-twist processes which require polarized quark and gluon distributions as 
input. Our input parameters have been adjusted to be compatible with global pa- 
rameters available current experiments. The values can be refined as further and 
more precise polarization experiments become available. A more precise parame- 
terization should also take into account corrections from QCD evolution, although 
this effect is relatively unimportant for helicity-dependent distributions. Our cen- 
tral observation is that the shape of the distributions is almost completely predicted 
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when one employs the constraints obtained from general QCD arguments at large 
X and small x. 

A remarkable prediction of our formalism is the very strong correlations be- 
tween the parent hadron helicity and each of its valence-quark, sea-quark, and 
gluon constituents at large light-cone momentum fraction x. Although the total 
quark helicity content of the proton is small, we predict a strong positive correla- 
tion of the proton's helicity with the helicity of its u quarks and gluon constituents. 
The model is also consistent with the assumption that the strange plus anti-strange 
quarks carry 3.5% of the proton's momentum and —10% of its helicity. We also 
note that completely independent predictions based on QCD sum rules also imply 
that the three-valence-quark light-cone distribution amplitude has a very strong 
positive correlation at large x when the -u— quark and proton helicities are paral- 
lei'"'. 
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FIGURE CAPTIONS 

Fig. 1: Model predictions for the non-perturbative polarized Aq{x) = q^{x) — 
q~{x) and unpolarized quark xq{x) = x{q'^{x) + q~{x)) distributions in the pro- 
ton. The polynomial forms satisfy sum rule and dynamical constraints. The leading 
Regge behavior at a; — » has intercept a = L12. By definition both quark and anti- 
quark contributions are included. Comparison with the MRS D'q parameterization 

[27] 

for the unpolarized quark distributions are also shown, (a): u{x) distributions. 
(b) d{x) distributions, (c) s{x) distributions. 

Fig. 2: Predictions for the non-perturbative polarized AG'(a;) =G^{x) — G^{x) 
and unpolarized gluon xG{x) = x{G^{x) + G^{x)) distributions in the proton. 
The polynomial forms satisfy sum rule and dynamical constraints. The leading 
Regge behavior at a; ^ has intercept ag = 1.12. Comparison with the MRS D'q 

[271 

parameterization for the unpolarized gluon distributions is also shown. 

Fig. 3(a): Model prediction for the polarized helicity structure function of the 
proton compared with experiment. Full line: sum of all flavors; dashed: only up 
quarks; dotted: only down quarks; dash-dotted: only strange quarks. We have 
multiplied our prediction with a PQCD correction factor 1 — (tts/Tr) = 0.92. The 
data are from SLAC-EMC (closed circles), EMC (closed squares), SMC (open 
squares) and SLAC-E143 (diamonds). Fig. 3(b): same as (a) but for the neutron. 

[21 

The data are from the SLAC E142 experiment . 

Fig. 4: Polarized helicity structure function of the deuteron. The data are from 
Ref. [1]. We have multiplied our prediction from the sum of proton and neutron 
contributions by a D— state depolarization factor 1 — {3/2)uj]j with ujjj = 0.058 
and the PQCD correction factor 1 — (as/Tf) = 0.92. 
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